Radiation from elementary sources in a uniaxial wire medium 
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We investigate the radiation properties of two types of elementary sources embedded in a uniaxial 
wire medium: a short dipole parallel to the wires and a lumped voltage source connected across a 
gap in a generic metallic wire. It is demonstrated that the radiation pattern of these elementary 
sources have quite anomalous and unusual properties. Specifically, the radiation pattern of a short 
vertical dipole resembles that of an isotropic radiator close to the effective plasma frequency of 
the wire medium, whereas the radiation from the lumped voltage generator is characterized by an 
infinite directivity and a non-diffractive far-field distribution. 

PACS numbers: 42.70.Qs, 78.20.Ci, 41.20.Jb 



I. INTRODUCTION 

Wire media, generically defined as structured materials 
formed by arrays of long metallic wires, ^^'^ is perhaps the 
class of metamaterials whose effective response is better 
understood. Particularly, during the last decade a vast 
body of theoretical methods and analytical tools have 
been developed which enable characterizing the effective 
electromagnetic response of wire-based materials in dif- 
ferent scenarios with great accuracy. '^^-'^^ However, a bit 
surprisingly, the problem of radiation by localized exter- 
nal sources embedded within wire media, has only been 
cursorily discussed in the literature. ^^"^^ 

This gap can be explained in part because of the pe- 
culiar electromagnetic response of wire media, which are 
typically characterized by strong spatial dispersion in the 
long wavelength limit, ^ and this property greatly compli- 
cates the analytical modeling. In simple terms, a medium 
is spatially dispersive if the polarization vector at some 
generic point in space depends not only on the macro- 
scopic electric field, but also on the gradient of the field, 
and possibly higher order derivatives.^^ 

The objective of this work is to characterize the radi- 
ation properties of elementary external localized sources 
placed within a wire medium using an effective medium 
approach. Specifically, we are interested in the follow- 
ing two scenarios: (i) a short vertical dipole is embed- 
ded in the wire medium (Fig. la) and (ii) an external 
lumped voltage source is connected across a gap in a 
generic metallic wire (Fig. lb) . As will be detailed in Sec- 
tion II, these sources are modeled in terms of the Dirac 
delta function (Dirac-5). It is, however, important to 
make clear at the outset that an effective medium de- 
scription of the radiation problem is only possible if the 
source is localized on a larger scale than the characteristic 
dimension of the metamaterial (e.g., the lattice constant 
a, see Fig. 1). Hence, the short vertical dipole consid- 
ered here should be understood as some external current 
spread over a region of space whose characteristic diam- 
eter in the xoy-plane is larger or equal than a, but much 
smaller than the wavelength. Similarly, even though for 
the purpose of illustration and discussion we say that in 



FIG. 1: (Color online) Uniaxial wire medium formed by a 
square lattice of metallic wires oriented along the z-direction. 
(a) Excitation based on a short vertical dipole embedded in 
the wire medium, (b) Excitation based on a discrete voltage 
source connected directly at the center of one of the wires. 



case (ii) the voltage source is connected across the gap of 
a single wire, it is more accurate to imagine such source 
as an array of voltage generators, distributed over a re- 
gion of space whose characteristic diameter is larger than 
a, being each voltage generator connected across a gap in 
a metallic wire lying within the mentioned region. With 
the exception of the immediate vicinity of these sources, 
the solution determined with our theory (based on the 
Dirac-5 distribution) should describe accurately the radi- 
ated fields. 

One of the challenges in the characterization of the ra- 
diation by a localized source within a wire medium is re- 
lated to the calculation of quantities such as the Poynting 
vector or the radiation intensity (i.e. the power radiated 
per unit of solid angle). Indeed, in general the usual form 
of the Poynting vector, S = E x H, does not hold in case 
of spatially dispersive materials. ^'^'^'^ Moreover, there is 
no known theory to determine the Poynting vector in a 
general spatially dispersive material, and the only case 
that is actually understood and for which closed analyt- 
ical formulas are available, is when the electromagnetic 
fields have a plane wave type-spatial variation. ^'^ In this 
work, we derive closed analytical formulas that enable 
calculating explicitly the Poynting vector and the elec- 



2 



tromagnctic energy density in uniaxial wire media for ar- 
bitrary electromagnetic field distributions. This is one of 
the key results of the paper. 

To do this, we rely on the theory of our earlier 
works/^'^"'' where wc have shown that the effective 
medium response of the wire medium can be modeled us- 
ing a quasi-static model, based on the introduction of two 
additional variables I and (p^ ■ What is remarkable about 
such a model, is that the relations between the macro- 
scopic electromagnetic fields and the additional variables 
I and are local in space. Therefore, such formalism 
enables describing the unconventional electrodynamics of 
the wire medium using a local approach, without requir- 
ing the definition of an effective spatially dispersive di- 
electric function, which would lead to nonlocal relations 
between the polarization vector and the macroscopic elec- 
tric field. ^ It is important to mention that the introduc- 
tion of the additional variables is not just a trick that 
simplifies the modeling of the wire medium: it is actually 
full of physical significance, and elucidates about the in- 
ternal physical processes that determine the macroscopic 
response of the metamaterial. Indeed, the variable I can 
be understood as the electric current that flows along 
the metallic wires (interpolated in such a manner that 
it becomes a continuous function defined in all space), 
whereas the variable can be understood as the aver- 
age potential drop from a given wire to the boundary of 
the respective unit cell (the potential is interpolated in 
the same manner as the current). For more details, the 
reader is referred to Refs. 17,25. 

This paper is organized as follows. In Section II, we 
briefly review the quasi-static model of the wire medium, 
and formulate the radiation problem for the two excita- 
tions of interest. In Section III we solve the pertinent 
radiation problem in the spectral domain. First, we dis- 
cuss the general case of a stratified (along z) structure, 
and after this we analyze in details the particular case of 
an unbounded uniform structure. In Section IV, we show 
that for an unbounded uniform structure the flelds radi- 
ated by the elementary external sources can be as well 
directly determined from the nonlocal dielectric function 
of the metamaterial. After this, in Section V we derive 
a general Poynting theorem that expresses the conserva- 
tion of energy in wire media, and in Section VI we use 
these results to obtain the asymptotic form of the Poynt- 
ing vector in the far-field, as well as the directive gain, 
directivity, and the power radiated by a short vertical 
dipole. The conclusions arc drawn in Section VII. In this 
work, we assume that in case of time harmonic regime 
the time variation is of the form e-'"*. 



II. LOCAL FORMULATION BASED ON THE 
INTRODUCTION OF ADDITIONAL VARIABLES 

In Refs. 17,25 it was shown that the internal physi- 
cal processes that determine the macroscopic response of 
a wide class of wire media, are intrinsically related to 
the dynamics of the electric current / along the wires 



and the additional potential whose physical mean- 
ing was already discussed in Introduction. In particu- 
lar, it was proven that for the case of straight wires ori- 
ented along the z-direction the macroscopic electromag- 
netic fields satisfy 

V X E = -jwMoH, (1) 

V X H = Je^t + + jw£/»E, (2) 

where eh is the permittivity of the host material, Ac = 
is the area of the unit cell, and a is the period of the 
wire medium (Fig. 1). Notice that unlike in our previ- 
ous works, "'^'''^^ here we admit the possibility of an exter- 
nal distributed current source iext- The electromagnetic 
fields are coupled to the current / and additional poten- 
tial ifw via a set of transmission line-type equations: 

^ = -{Z^+ jcoL)I + E, + Vc^tAcS {x, y, z) , (3) 

^ = -jujC'(pw (4) 
oz 

In the above, C, L, and Zw represent the capacitance, in- 
ductance, and self-impedance of a wire per unit of length, 
respectively, and explicit formulas for those parameters 
can be found in our previous papers. The real part of 
is related to the ohmic loss in the nicitallic wires, whereas 
its imaginary part is related to the kinetic inductance of 
the electrons in the metal. As compared to Refs. 17,25 
now we allow for an external lumped voltage source (with 
amplitude Vext) to be placed across a gap in the wire in 
the central unit cell. It is simple to check based on the 
theory of Ref. 17, that this lumped voltage source is mod- 
eled by the term Vext^cS {x, y, z). Notice that similar to 
the current and additional potential, the lumped genera- 
tor is interpolated so that it becomes a function defined 
over all space. 

In the next section, we will determine the solution of 
the radiation problems sketched in Fig. 1, based on the 
system of Eqs. (1)~(4). It should be emphasized that such 
a system of equations is local in the sense that all the rel- 
evant medium parameters (C, L, Z^, Sh, and /xq) are 
independent of the gradient and higher order derivatives 
of the electromagnetic fields, /, and ip,^]. This contrasts 
with the usual formulation based on the effective dielec- 
tric function, which does not require the introduction of 
additional variables hut in which the dielectric function 
depends explicitly on the wave vector. This will be 
further discussed in Section IV. 

For future reference, we note that from Eqs. (3)-(4) it 
follows that 

= -juj [E, + VextAcS (x, y, z)] . (5) 

In the above, it was supposed that L, C, and Zyj may 
depend on z (but not on x and y), which can happen in 
case of a stratified wire medium (with direction of strati- 
fication along z), such that either the permittivity of the 
host medium or the radii of the wires varies with z. 
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III. THE RADIATION PROBLEM 

Next, we derive the solution of the radiation problem 
in terms of two scalar potentials. Wc admit that the 
external current density describes a short vertical dipole, 
so that that Jext = ji^PeS {x, y, z) z, where Pe represents 
the electric dipole moment. Since the Maxwell equations 
are linear it is possible to solve the two radiation problems 
sketched in Fig. 1 simultaneously. This will be done in 
what follows. 



applicable t o both square and hexa gonal wire lattices 
kp ~ (l/a)-y/27r/ log[a2/4r^(a — r^,)], where is the ra- 
dius of the metallic wires. Using the fact that for straight 
unloaded wires LC = i^o^h, it follows that Eqs. (8) and 
(9) are equivalent to 

d 1 9$ 

"a^ + ^* * + ''^^ + ''^^ " '^"^ ^' ' ^^^^ 



(13) 



A. Solution in terms of two scalar potentials for 
the general case of a stratified structure 

For generality, in this subsection, we admit that L, C, 
Zy,, and eh may depend on z, which as discussed previ- 
ously, may be useful to study problems of radiation in 
stratified media. We look for a solution of (l)-(4) such 
that the macroscopic electromagnetic fields are written 
in terms of a Hertzian potential $ so that 



H = V X {jw^z} , 



(6) 
(7) 



It can be easily verified that (6)-(7) satisfy, indeed, the 
Maxwell equations (l)-(2), provided that 



d 1 (9$ 

^'^^ -^+^t^ 

OZ Sh oz 



= -peS {x,y,z) , 



(8) 



where wc put kf^ = w^/ioe/i and = —jtoZyjC. Hence, to 
determine the solution of our problem, wc need to solve 
this coupled system of partial differential equations with 
unknowns $ and V'- 

To do this, it is most convenient to work in the Fourier 
domain. Defining $ and ip as the Fourier transform of $ 
and V' in the xy plane, respectively, so that 



^ = j j ^e^''^-'=+^yy^dxdy, 



(14) 



and ij) is defined similarly, it follows that 
rl 1 (9<t> 

e/.^-^ + (fc^ - kl) $ + fc> = -Pe5 {z) , (15) 

^^-zhb' + + + kf^ = PefS (z) , (16) 

where k^ = k'^ + ky. Thus, we have reduced the radiation 
problem to the solution of a system of linear ordinary 
differential equations. 



o2 o2 

^ + Hence, substituting (7) into (5) 



where 

and using the above result, it follows that 



jujAc 



C 

AcEh 

C 

AcSh 



[ShE^: + EhVextAcS {x, y, z)] 

-Vt^-PefS{x,y,z) 



(9) 



where we defined the effective dipole moment for the com- 
bined excitations: 



Pef =Pe- ShAcVext ■ 



(10) 



Notice that pe/ depends on both pe and Vext , because we 
allow for the simultaneous excitation of the wire medium 
with the two pertinent types of elementary sources. 

For convenience, let us introduce the auxiliary poten- 
tial 



1 



k^ jwAc ' 



(11) 



where kp = fi^/ {LAc) is the so-called plasma wave 
number of the wire medium,^'^^'^'' which may be cal- 
culated using, for example, the approximate formula 



B. The case of a homogeneous medium 

From hereafter, wc restrict our attention to the par- 
ticular case of a homogeneous and uniform medium, for 
which the structural parameters Sh, C, L, and Z^, can 
be assumed independent of z. In such a case, the system 
(15)-(16) can be rewritten in a compact matrix notation 
as follows: 



dz^ \ ip 



A-2 



kp 

_h.2 1,2 _ iL2 _ o2 
'^t '^p l^h Pc 



-Pe 
Pef 



(17) 



The general solution of the homogeneous problem, when 
Pe = Pef = 0, can be easily found using standard meth- 
ods, and is given by 



= (C+e-T^'^^-FCfe+''^ 



Mz\ \ kp 

ll - iTM 



(18) 



where 7^ = /e^ — /e^, Ci^ with i = 1,2 are integration con- 
stants, and 7gT and 7tm are the propagation constants 
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along the z direct ion of the so-cahed quasi-transverse elec- 
tromagnetic (gT) and transverse magnetic (TM) modes 
supported by the bulk wire medium. These parameters 
are defined consistently with Refs. 12,26, and satisfy 



By direct substitution into (17), it is readily found that 
the unknown constants CqT and Ctm are required to 
satisfy 



1 



7TM =j kl-- {kl + kf - (3^ + 

7,T=J kl-\{kl+kl-Pl- 
^(kj + k^-P^y+AkfPl 



(19) 



(20) 



This yields 



kp 



IqrCqT \ 2 2 



^2 

+ ITmCtM ( 2^2 

Ih ~ ITM 



1 f -Pe 

2 V Pef 



(22) 



In case of perfectly conducting wires, we have Zy^, = 0, 
and thus Pc = 0. In such a case the propagation constants 
of the qT and TM modes reduce to the well known forms, 

fc^, respectively.^ 



jU2 

ftp 



k? 



IqT = jkh and 7tm = 

Since the solution of (17) is obviously an even function 
of 2;, we may try a solution of the form 



i^J ' \^ h- iIt 



+ Ctms ■ \ 2 '2 



k?: 



7rM|z| I '"p 



(21) 



CqT = 



Ctm = 



1 hh-^TM)Pe + klPef 1 

279T IqT - Itm kl ' 

1 {lh-llT)Pe + klpef 1 



27TM 



7tm 7gT 



kl' 



(23) 
(24) 



Substituting this result into Eq. (21), we finally obtain 
the desired solution 



'0 / 27gT 



1 hi - Itm) Pe + kip, 



7gT Itm 



1 

— 1 



1 {ll - iIt) Pe + fcpP. 



27TM Itm - %t 
I 



1 

"fcT 



e-^""l^l 7g-7^M 1 • (25) 



The inverse Fourier transform of $ is given by 



C. Perfectly electric conducting wires 



$ = 



+ 00 

^ j ^Jo{ktp)ktdkt 



(26) 



where Jq is the zero-order Bessel function of the first kind, 
p = \/ + 2/^) and in the second identity we used the fact 
that $ is a function of kt- In general, this Sommerfeld- 
type integral can only be evaluated using numerical meth- 
ods. Obviously, it is possible to write a similar formula 
for ip. 



Let us now study what happens when to a first approx- 
imation the metal can be modeled as a perfect electric 
conductor (PEC), so that w 0. In such a situation, 
Eq. (25) simplifies to 



1 j^2 Pe - Pef ^--.^M 



2lqT "k^ + k^ 
1 kfpe + k^Pef 



2^TM 



kl + fc? 



-1tm\A 



(27) 



We will discuss separately the two scenarios of interest. 
Let us consider first that Vf-^t = 0, so that the metama- 
terial is excited solely with the short vertical dipole. In 
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this case, Pef = Pe, and thus we obtain simply 

(*)^-247=-™'='(-i)- P«) 

The corresponding inverse Fourier transforms can be eval- 
uated analytically in a trivial manner. This yields 

with kef = kf^ — fc^ and r — + y'^ + z^. 

Hence, it is immediately seen that the qT mode (which 
in the case of a lossless metal is exactly transverse electro- 
magnetic mode [TEM] with respect to the z-direction^) 
does not contribute to the radiation field of the short ver- 
tical dipole. This may look surprising at first, but it is 
actually simple to understand. Indeed, it is well known 
that the electric Green dyadic in a periodic structure (e.g. 
a photonic crystal or a metamaterial) can be written as a 
weighted summation of terms such as E„ (g) E„ where E„ 
stands for a generic natural mode of the system, and ® 
represents the tensor product. In particular, this im- 
mediately implies that a TEM mode (with respect to the 
z-direction) cannot possibly contribute to the field radi- 
ated by a short vertical dipole, because its contribution 
would be proportional to Etbm (Etba/ • z), whereas for 
a TEM mode Et_bm -2 = 0. Notice that this discussion 
applies actually to the microscopic electromagnetic fields 
(before homogenization on the scale of the lattice con- 
stant), but it clearly indicates that the TEM mode cannot 
contribute as well to the radiated field in the framework 
of a macroscopic theory, consistent with Eq. (29). It is 
interesting to note that in the presence of loss the con- 
tribution of the qT mode to the radiation field does not 
vanish [the first addend of Eq. (25) does not vanish when 
Pe = Pef], which is fully consistent with the microscopic 
theory, because in case of loss the electric field associated 
with the qT mode has a small longitudinal component 
(i.e. a component along the z direction). 

The result (29) implies two unexpected things. First, 
despite the anisotropy of the wire medium, the wavefronts 
are spherical surfaces! Second, the emission of radiation 
is possible only above the effective plasma frequency of 
the metamaterial, LOp — kp/ ^fshjio. The latter property 
is actually a direct implication of the TEM mode not 
being excited, as discussed above. To illustrate variation 
in space of the potential $, we plot in Fig. 2 the contour 
plots of $ for two different frequencies of operation. The 
wire medium is formed by PEC wires with r^^ = O.Ola 
standing in a vacuum. The plasma wave number of the 
effective medium is kp — 1.38/a. Thus, the example of 
Fig. 2a corresponds to a frequency below LOp, whereas 
the example of Fig. 2b corresponds to a frequency above 
LOp. This explains that in the former case the potential 
is strongly localized in the vicinity of the dipole, whereas 
in the latter case the potential decays much more slowly 
as 1/r. 

At first glance, the result (29) could suggest that the 
electric field radiation pattern should be similar to that 
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FIG. 2: (Color online) Contour plots of the amplitude of the 
potential <^ (arbitrary logarithmic unities), (a) uia/c = 0.5. 
(b) uia/c — 1.5. The wire medium is formed by PEC wires 
with = O.Ola standing in a vacuum, and is excited by a 
short vertical dipole. 

of a Hertzian dipole standing in a homogeneous isotropic 
plasma with a Drude dispersion. As it will be shown in 
Section VI, this is not true. 

Next, we consider the case Pe = 0, so that the wire 
medium is excited by a lumped voltage source (Fig. lb). 
In this scenario, the contribution of the TEM mode to 
the radiation field does not vanish. Indeed, the inverse 
Fourier transform of the first term in the right-hand side 
of (27) can be readily calculated and is equal to 

%T = -Pe/;r^e-^«-l^l^ifo (kpp) , (30) 
ZjqT 27r 

where Kq is the modified Bessel function of the second 
kind. On the other hand, the auxiliary potential (11) as- 
sociated with the wire current satisfies V'gT = —-^^t^qT- 

The result (30) is quite remarkable, because it predicts 
that the Hertzian potential, and hence the electromag- 
netic fields, varies with z simply as e"'''''^'^! = e^-''""'''^! , 
and hence the radiated field is simply guided along z, 
without any form of decay. Moreover, is strongly 
localized in the vicinity of the z-axis, within a spatial 
region whose characteristic diameter is determined by 
Ap = 27r/fcp. It should be mentioned that $gT is actually 
singular over the z-axis (it has a logarithmic singular- 
ity). Such a singularity occurs because of the adopted 
(S-function model for the lumped voltage generator (also 
because the macroscopic model effectively assumes in- 
finitesimal wire separation: o — > 0). The singularity dis- 
appears if one considers a less localized model for the 
discrete source, e.g. if S (x, y, z) is replaced by g [p) 5 (z), 
where g is some function of p concentrated near the ori- 
gin. Even for such a source, the electromagnetic fields 
are characterized by a non-diffractive pattern. This is 
explained by the "canalization" properties of the wire 
medium, which enable the transport of the near-field with 
no diffraction. As far as we could check, the inverse 
Fourier transform of the second addend of (27), i.e. the 
contribution of the TM mode when Pf, = 0, cannot be 
written in terms of the standard special functions, and 
hence it needs to be calculated numerically using (26). 
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particular case of the uniaxial wire medium formed by 
straight wires, the effective dielectric function is^'^^-^'^ 
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FIG. 3: (Color online) Contour plots of the amplitude of the 
potential $ (arbitrary logarithmic unities), (a) a;a/c = 0.5. 
(b) bjajc = 1.5. The wire medium is formed by PEC wires 
with = O.Ola standing in a vacuum, and is excited by a 
lumped voltage generator. In case (b) an interference pattern 
of the TEM and TM modes is observed. 



In Fig. 3, we plot the contour plots of $, for the same 
example as in Fig. 2. It is seen that when wa/c — 0.5 
(Fig. 3a), i.e. below the effective plasma frequency, the 
fields are strongly concentrated close to the z-axis, and 
are guided away from the source with no diffraction. The 
contribution from the TM mode appears to be residual. 
On the other hand, above the plasma frequency (Fig. 3a), 
there are clearly two distinct emission channels, one as- 
sociated with the TEM mode and another with the TM 
mode. 



IV. NONLOCAL DIELECTRIC FUNCTION 
APPROACH 

The objective of this section is to prove that the ra- 
diation problem in an unbounded uniform structure can 
be as well solved using the standard nonlocal dielectric 
function formalism. ^^2,17 rpj^g ^^^^ ^ stratified struc- 
ture, which can be easily handled with the theory of Sec- 
tion III, is out of reach of the nonlocal framework (it 
could however be handled with a combination of mode 
matching and additional boundary conditions, "'^"'^'^^ but a 
detailed discussion of it is out of the scope of this paper). 

In a spatially dispersive medium, the Maxwell equa- 
tions may be written in a compact form in the space 
domain as follows: 

V X E = -jcj/ioH, (31) 

V X H - Je,t + jcjf (w, jV) • E, (32) 

where the dyadic operator e(a;, jV) represents the effec- 
tive dielectric function of the material. Notice that in the 
space domain the effective dielectric function should be 
regarded as a function of the gradient V. This contrasts 
with the formulation of Section II where all the structural 
parameters are independent of V. It is also possible to 
write the term e (w, j V) • E as a spatial convolution.^'^ In 
the spectral (Fourier) domain, in which j V k, in the 



(33) 



where ^ — [2^,/ L),JehJlQ, and the rest of the symbols 
are defined as in Section II. Notice that the effective 
dielectric function depends explicitly on ^ J ^ • 

Despite the apparently complicated form of Eqs. (31)- 
(32), the radiation problem can be readily solved in the 
spectral domain in case of an unbounded uniform struc- 
ture. Indeed, by calculating the Fourier transform of both 
sides of the equations (31)-(32) with respect to all the 
space coordinates, so that jV O k, it is readily found 
that 

k X E(a;,k) = w//oH(w,k), (34) 
k X H(a;, k) = -wf (w, k) • E(a;, k) - ujV^^t{uj, k), (35) 

where jwPe^t (w, k) — Jea;t(a;,k) is the Fourier-trans- 
formed source term. After some straightforward manipu- 
lations, we find that the Fourier transform of the electric 
field is 



E = jw/io w^^oE (w, k) -(- kk - k^l 



•J. 



(36) 



and hence the electric field in the space domain can be 
formally written as 

E(r)^^x 
(2^)' 

wVo? (w, k) + kk - fc^fj • Je^t (k) e^^'^-'^d^k. 

(37) 

Notice that, at least a priori, in the nonlocal dielec- 
tric function framework we can only consider excitations 
based on an external density of current (Fig. la). The 
characterization of the excitation based on a lumped volt- 
age source requires the knowledge of internal degrees of 
freedom of the wire medium (e.g., the current along the 
wires and the additional potential), which are not de- 
scribed by the effective medium model. Nevertheless, 
ahead we will show that a lumped source Vext can also 
be modeled by a suitable equivalent 3ext- 

Next, we obtain the solution of the radiation prob- 
lem when Jea;t (r) = jujpez6{r), or equivalently when 
Pej;t(a;,k) — peZ,. Instead of attempting to calculate the 
integral (37) directly, we will solve instead (34)-(35) by 
introducing the Hertz potential He. In this manner, we 
write the Fourier-transformed fields as follows: 



E 
H 



W^E/i/ZoIIe - k(k • He), 
UEhk X lie. 



(38) 
(39) 



This form immediately satisfies (34). From (35) we find 

k X (k X He) = -— • [klU, - k(k • He)] - (40) 

Sh £h 



7 



Using (33), after some trivial vector algebra, we obtain 

{kl - e)n, = -x^^zz ■ [kin, - k(k • He)] - (41) 



where 



Hp 



kl - j^kh - kl ' 



(42) 



Calculating the vector product of (41) by z we find that 
z X He = 0, (43) 
and therefore, n,. = II^z. From this and (41), 



[kl-e + xzz{kl-kl)] n, = - 



eh' 



For PEC wires ^ = 0, and thus from (42) we have 



eh 



from which 



n^(a;,k) = —2 — 12 — 

ehyi^fi lip n ) 



and, thus. 



n,(u;,r) 



Pe e 



Sh 47rr 



(44) 



(45) 



(46) 



(47) 



which is the same as (29), because $ = e/ill^. 

In the general case in which the metal has a plasmonic- 
type response, ^ ^ 0. Introducing the notation = 
—j^kh we obtain from (44) 



n^(w,k) 

Pe 



Eh {kl - k^){kl + - fc2) - fc2(fc2 _ ^2) 



Pe kl + _ fc2 

eh {kl + lTM){kl + iItY 



(48) 



where 7tm and are given by (19) and (20). 

Calculating the inverse Fourier transform with respect 
to kz we find 



I\.z{uJ,)s.f,z) 



Pe ( K+PI+^It ^ 

2eh [jqH'y'^T - I'tm)' 



-7<,t|z| 



_|_ K+ Pc+ ITM ^—iTM\z\ 

1tm{1tm - Igr) 



(49) 



At first glance this result looks different from (25), but 
one may verify that kl + + 7^^ = (7^ - 72^) + k^, 

and, similarly, kl + /S^ + = {ll " iIt) + ^p- Thus, 
we recover (25) with = p^f, which corresponds to the 



case Vext = 0, consistent with our assumptions in the 
beginning of this section. 

Surprisingly, the lumped voltage source Vext in (3)- 
(4) can be equivalently represented within the nonlocal 
dielectric function model with some distributed current 
density 3ext,v in the unbounded wire medium. To show 
this, we consider the Fourier-transformed equations (3)- 
(4) (for simplicity we let =0), from which the Fourier- 
transformed current I{oj, k) can be expressed as 



/(w,k) = -jojShA, 



° i.2 _ 1.2 



[Ez{LU,k) + VextAc]. (50) 



When this expression is substituted into the Fourier- 
transformed equations (l)-(2), the i?z -proportional term 
of (50) is combined with the term joJEhE which results in 
the spatially dispersive permittivity (33), and the Vext- 
proportional term occurs as an additional external cur- 
rent density 

k'^A 

Jext,y(t^,k) = -jwShVext ,2^ L ^" (^^) 

l^h '^z 

Therefore, applying the inverse Fourier transform, we find 
that 

klA^Vexte-^^'^^''^ 
Jea:t,v'(w,r) = -^^ 5{x,y)i, (52) 



where r]h = \/ Ho/ en- Thus, a lumped voltage source in- 
serted into a wire of the unbounded uniaxial wire medium 
(with PEC wires, = 0) may be equivalently repre- 
sented with a line of z-dircctcd wave-like current (52). It 
is curious to note that while the lumped voltage source 
excitation is localized at the origin, the equivalent cur- 
rent density is distributed over the entire z-axis. At first 
sight, this may look inconsistent with causality. However, 
it is simple to verify that such a current is just a wave 
emerging from the discontinuity point at z = 0. Indeed, if 
one calculates the inverse Fourier transform of (52) with 
respect to time, it is found that: 



A ~ 

^ext,V {t,r) = kl—^Vext 
^Vh 



t-^]zS{x,y), (53) 

Vh/ 



with Vh = If-s/ehfJ-o the velocity of propagation in the 

host material, and Vext {t) the inverse Fourier transform 
of Vext (<^)- The above formula is manifestly consistent 
with causality, because the excitation at a given point z 
only depends on the excitation at the origin with a delay 
\z\/vh. 



V. ENERGY CONSERVATION IN THE 
UNIAXIAL WIRE MEDIUM AND POYNTING 
THEOREM 

In what follows, we prove that the formulation based 
on the introduction of additional variables of Section II, 
enables formulating an energy conservation theorem and 
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the definition of a Poynting vector in the uniaxial wire 
medium. 

Wc start with equations (l)-(4) written in time do- 
main. The host permittivity Sh is assumed dispersion- 
less and lossless, and the wires are modeled by a self- 
impedance of the form Z.u;{uj) = jujLi^in + R, where the 
parameters Lkin and R are independent of frequency. For 
metallic wires with radius standing in air and de- 
scribed by the Drude model with the plasma frequency 
ujm and the collision frequency F, these parameters are 
Lkin = l/(eo7rr^w^) and R = r/(£o7rr^w^). 

Thus, in the time domain the equations (l)-(2) and 
(3)-(4) may be written as 



V X E = -fxo 



dU 



V X H = £^ 



dt ' 



' ext-> 



(54) 

(55) 



and 



— --(i + iki„)^ 
dl_ _ ^dip^ 

dz at ' 



RI + E,+£ext, (56) 
(57) 



where Sext is the effective EMF of the voltage sources 
inserted into the wires, per unit length of the wires [e.g., 
for a lumped source V^xt inserted into a wire at r = 0, 

Sext = VextAcS{r)]. 

Following a standard procedure, we obtain from (54)- 
(55): 



V.[ExH] = -| 



-E-Je,t-^. (58) 



and Pcxt as the volume density of the power transferred 
by the external sources to the medium. In the absence 
of loss, i.e. when i? = 0, the term W is univocally iden- 
tified with the density of stored energy. In contrast, if 
loss is present, then it is generally impossible to separate 
the energy storage rate from the energy loss rate when a 
metamaterial is considered macro scopically. 

However, if the micro structure of a metamaterial is 
known, the stored energy can be found from a consistent 
physical model that fully describes the processes within 
a unit volume of the metamaterial. Thus, if we assume 
that the Drude model is such a consistent model for the 
dynamics of the free electron plasma in metals, then (62) 
preserves the meaning of the stored energy density even 
when i? > 0. In this case, the quantity Pioss has the 
physical meaning of an instantaneous power loss density. 

Evidently, in time-harmonic regime the time-averaged 
Poynting vector is given by 



S„. = iRejExH* + ^z 



(65) 



It can be checked that in the lossless case (Re {Z^,} = 0), 
and for the case of fields with a spatial dependence of 
the form e~^^'^ with k real-valued this reduces to the 
formula, 

Sa.,/ = iRe{(E x H*)J - I^* • ^ • ^' (66) 

with / = .T, y, z and e defined as in Eq. (33), which is ap- 
plicable to plane waves in general lossless spatially disper- 
sive media.^^'^^'^^ The application of the above formula 
to wire media has been considered in several works. 



On the other hand, from (56)-(57) we have 



dz 



d_ 
'dt 



- Ri' + Ej + e^xti, 

(59) 

where itot = L-\- ikin- 

Diving the last relation by and adding it to (58) we 
obtain the conservation law 



_ dW 

V • S = -Ploss + -Pext) 



where 



S = E X H 
W 



2 
RP 



2Ar, 



2Ar. ' 



-E J 



ext ■ 



(60) 

(61) 
(62) 
(63) 
(64) 



The vectorial quantity S in (60) and (61) may be under- 
stood as the Poynting vector in the uniaxial wire medium, 



VI. RADIATION PATTERN IN THE PEC CASE 

Next, we obtain the radiation pattern, directive gain, 
directivity, and radiation resistance for the case of a short 
vertical dipole radiating in a wire medium formed by PEC 
wires (Fig. la). We do not discuss in details the case 
wherein the metamaterial is excited by a lumped voltage 
source, because as discussed in Section IIIC in such a 
scenario the radiated field is guided along the z-axis with 
no decay. In particular, this implies immediately that the 
directivity in such a configuration is infinite. 



A. Asymptotic form of the radiated fields 

To begin with, we obtain the asymptotic form of the 
field radiated by a short vertical dipole iyext = 0) embed- 
ded in a wire medium formed by PEC wires when r — > oo. 
Evidently, from the results of Section III C, unless the fre- 
quency of operation is larger than the plasma frequency 
of the effective medium the radiated fields will decay ex- 
ponentially away from the source. Hence, in what follows 
we assume that w > = kp/ ^/Jloeh, so that k^f > in 
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Eq. (29). Substituting (29) into Eqs. (6)-(7), it can be 
easily checked that 



H= — ujkef^ sinO If 



(67) 



£h 



1 



sin 06* 



indicates that the identities are 
$ = PpT^.e~^^''f'', and (^r,9,ip^ 



where the symbol 

asymptotic (r ooj, <P = pe^^, 

define an orthogonal reference system associated with the 
usual spherical coordinate system {r,9,ip). As seen, un- 
like what happens in an isotropic medium, the electric 
far field has a radial component. The amplitude of the 
electromagnetic fields varies asymptotically as l/r, and 
Ee = T]efH^ with rjef = ujfj^o/kef. 

Similarly, substituting (29) into Eqs. (3)-(4) and (11), 
it is found that the asymptotic forms of the current and 
additional potential are 



= /— ^: COS( 



jk, 



COS0$. 



(69) 
(70) 



Thus, from (68) and (70), we see that the time averaged 
Poynting vector (65) in the far-field is 



1 k^ ^ ^ 

Sat; = T^w^Mo I'i'l^ ^e/ sin 6' -f COS 6* ^ + sin6'f 

+ i^fc>cos6'|$|^z. 

^ £h 



(71) 



Straightforward calculations show that the Poynting vec- 
tor only has a radial component: 



1 / A: \ 

Sav = ^kefio^HQ |$|' f sin2 6 + -^ cos2 j f. 



(72) 



Hence, in part surprisingly, it is seen that a short vertical 
dipole embedded in a wire medium can radiate energy 
along the direction of vibration, i.e. along the z-axis! 
Moreover, in the limit cj Wp the radiation pattern be- 
comes isotropic: Sat, ~ ifce/a;"^/io |<i>|^ f . Note, however. 



that for oj 



We 



we have fee/ = and thus the Poynting 



vector vanishes in the far-field. However, slightly above 
Up the emission of radiation is certainly possible. Notice 
also that in the limit where /cp — >■ 0, we recover the far- 
field of a short vertical dipole embedded in a dielectric 
with permittivity Eh- 



B. The radiation intensity, directive gain, 
directivity, and radiation resistance 



The radiation intensity of the short vertical dipole, U 
limr^oo r'^Sav, is given by 



U = ^fce/wVo f sin^ e + ^ cos^ 6 



(73) 



Hence, the power radiated by the dipole, Pr, 
2tt J U sin 6d9, is such that 



/ udn 



Pr. 



(74) 



The directive gain, g = AnU/Prad, is 
9(0,^) 



(75) 



Since kh > kp it can be checked that the direction of 
maximal radiation is — tt/2. The directivity of the 
short vertical dipole is, thus. 



D 



2+k^p/kr 



(76) 



which therefore increases from unity (for kh ~ kp) up to 
3/2 in the limit kh 3> kp. 





Z ^ 






— 






A;,, =bmX, 




^^^^ '^N^^N 


'III \ 
ill \ 




; \ y \ \k, = lo.ok / 

: ^^^^^^ 


^ ^'/,=2-5o^„ ) J ; 
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FIG. 4: (Color online) Polar plot of the directive gain of a 
short vertical dipole embedded in the uniaxial wire medium 
for different frequencies of operation (lj = kh / y/shjlo) ■ 

In Fig. 4 we show a polar plot of the directive gain of 
the short vertical dipole for different frequencies of oper- 
ation, normalized to the effective plasma frequency. In 
agreement with the previous discussion, it can be seen 
that the radiation pattern becomes more directive for in- 
creasing values of the frequency, and that for lu k ujp the 
radiator resembles an isotropic radiator. 

To conclude, we note that if the dipole is fed by a 
current /q and has infinitesimal height c?Z, then the cor- 
responding dipole moment is such that wjpel — 
Thus, it follows that the radiation resistance {Rrad ~ 
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'^Pradl l-^ol ) of such an elementary source is given by 

= ..(§^.e,fc.(l + ^), (77) 

where rjh = ^/jJ-o/^h is the impedance of the host mate- 
rial. 



VII. CONCLUSION 

In this work we have studied the radiation of two 
types of elementary sources embedded in a uniaxial wire 



medium and derived a general energy conservation theo- 
rem. The main challenge of the radiation problem is re- 
lated to the mctaniatcrial being spatially dispersive. We 
have shown that the radiation problem can be solved by 
considering either a nonlocal dielectric function frame- 
work or, alternatively, a local model framework based on 
the introduction of additional variables. However, only 
the latter approach enables considering stratified media 
and calculating quantities such as the Poynting vector or 
the directive gain. It was shown that the emission of ra- 
diation by a short dipole in a wire medium has several 
anomalous features, such as a uniform directive gain near 
the effective plasma frequency. On the other hand, the 
radiation by a lumped voltage generator results in a non- 
diffractive beam that is localized in the vicinity of the 
^;-axis, and corresponds to an infinite directivity. 
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